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Joule -Thomson Process in the Liquefaction of Helium 1 

Edmund H. Brown and John W. Dean 

A dimensionless, normalized, correlating function is introduced for the specific enthalpy 
of helium. Using this function, the consistency of various enthalpy data is determined. 
These data are then used to obtain curves of helium liquefaction yield in the region of variable 
specific heat in terms of heat-exchanger efficiency. 



1. Introduction 

Many different devices and cycle arrangements 
can be used in partially cooling gas for helium lique- 
faction. With the exception of the Simon free- 
expansion type liquefier, however — which is suitable 
only for small-yield, batch production — all such 
arrangements have a final stage at the lowest tem- 
peratures, in which a so-called Joule-Thomson process 
takes place. This process may be considered to in- 
clude cooling of the high-pressure gas by the un- 
liquefied return gas in a "Joule-Thomson heal 
exchanger'', as well as cooling of the high-pressure 
gas by the Joule-Thomson effect, that is, isenthalpic 
expansion through a "Joule-Thomson valve." If 
the initial state of the high-pressure gas at the top 
of the J-T heat exchanger (which is most con- 
veniently defined by its pressure and temperature) 
is fixed within a certain region, a definite mass frac- 
tion a of the gas passing through the J-T valve will 
be liquefied. 

To determine the yield a of a liquefier, the specific 
enthalpies at certain points in the J-T process must 
be known as functions of P and T. The most readily 
available data on the specific enthalpy, h, of helium 
is from the monograph of Keesom [l], 2 or the tem- 
perature-entropy diagram of Keesom published by 
the International Institute of Refrigeration [2]. 
The errors and inconsistencies of these sources below 
about 20° K, however, have been pointed out in a 
series of papers by Zelmanov [3, 4, 5, 6]. In addition, 
Zelmanov [6] has recalculated and redrawn both 
T-s and h-P diagrams which, as indicated later, 
appear to be of good accuracy and consistency. 
The Zelmanov T-s diagram is reproduced by 
Zemansky [7]; and the h-P diagram, by Collins [8]; 
who, however, also reproduces the Keesom T-s 
diagram . Because of the crowded nature of the iso- 
bars in the region of interest for the J-T process on 
a T-s diagram, the h-P diagram is more convenient 
for liquefier calculations. 

At the beginning of this study there was no easy 
way to judge the relative value of various conflicting 
sets of data. In order to determine, at least, the 
consistency of such data, a new dimensionless, 
normalized correlating function w(P, T) was intro- 
duced. With this function, the consistency of 
Zelmanov's work was evident, and both the lack of 
self-consistency of the T-s diagrams of Keesom, and 
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the inconsistency of these diagrams in certain regions 
with Keesom 's own density data was apparent. 
In addition, it appears that there may be some value 
in this new function w(P, T) in theoretical studies, 
since it implies that, above certain pressures, the 
enthalpy of a real gas may be approximately given 
by the sum of two linearly independent functions, 
one of the pressure alone, the other, of the temper- 
ature alone. However, not enough time was avail- 
able in this study for further investigations in this 
direction. 

A second factor which complicates the design of 
helium liquefiers is the behavior of the specific heat 
c p of the real gas at low temperatures. At high 
temperatures c p of a gas gradually increases with 
increasing pressure; but, at a lower temperature, 
dependent on the pressure, this trend is reversed. 
Ignorance of this effect has occasionally led the 
unwary to design liquefiers in which the temperature 
difference AT at the bottom of the J-T heat ex- 
changer would have to be reversed. This phe- 
nomenon, which has been hinted at by Keesom and 
discussed qualitatively by Zelmanov, is investigated 
further in this study. In order to determine per- 
formance, a method of analyzing J-T heat exchangers 
in terms of one parameter — the heat exchanger 
efficiency rj — is presented, and graphs of yield curves 
for helium based on this method are included. 

The purpose of this paper, then, consists of three 
things: First, and most important, to affirm the work 
of Zelmanov (an adaptation of Zelmanov's h-P 
diagram is included as figure 1); second, to present 
the correlating function w(P, T) and discuss some of 
its properties; and third, to show quantitatively the 
effect of the variation of c p in the J-T heat-exchanger 
on the liquefied fraction a. 

2. Correlating Function w{P,T) 

The enthalpy of a gas can be expressed as 

h(P,T) = fc p (P,T)dT+G(P) 

where is a constant of integration. For an ideal 
gas, c p (0,T) = c° p ] and, for helium, which is a mona- 
tomic gas, c J p = (5R/2) = 5.193 j/g-° K. In defining 
the enthalpy of the ideal gas, h°=h(0,T), the con- 
stant of integration 6?(o) was equated with the latent 
heat of evaporation at 0° K, L°, which is equivalent 
to fixing the enthalpy (which at this point is equal 
to the internal energy) of the liquid at P=0, T=0 
as zero. According to Swenson [9], the extrapolated 
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value of L° is 14.95 j/g. Thus, the expression used 
for helium in the ideal gas state is 

A°=c;r+L°=5.193r+ 14.95 j/g. 

Just as an enthalpy in an "ideal gas" state can be 
defined by extrapolating to zero pressure, so an 
enthalpy in an "ideal condensed phase" state, 
ho=h(P,Q), can be defined by extrapolating to zero 
temperature. If V(P,0) is the specific volume at 
zero temperature, then 

(W T - =V(P > 0)=Vme 

An attempt was made to express the extrapolated 
isothermal compressibility at 0° K, k(P,0) = — (1/V)- 
[(Z>V)/(dP)] T , by an empirical equation which would 
have led to an expression for h as a function of ex- 
ponential integrals, Ei(k). However, since this 
required some simplifying assumptions, was not par- 
ticularly suited for engineering calculations, and 
lacked theoretical justification, h was finally de- 
termined by extrapolating Keesom's density data 
[1] with the aid of expansivity and compressibility 
data of Atkins and Edwards [10]. The resulting 
curves for ho, V(P,0), and also k(P,0), are included 
as figures 2, 3, and 4. It should be mentioned that 
Swenson's extrapolated value [9] of F(0,0), for the 
liquid seems to be in error. 

Having, now, both h°=h(0,T) and h =h(P,0), the 
remaining "anomalies" of both the real gas and the 
condensed phases can all be contained in a function 
w(P, T) defined by 

_ h(P,T)-ho(P) 

which is not only dimensionless, as is evident from 
the definition, but has a total range of zero to one; 
zero, for the condensed phases as T-^0° K; one, for 
the gas as T-^oo ; one, for the gas as P— >o; and zero, 
for the condensed phase as P— ><» . 

Here, we use the symbols T-^oo or P-^oo merely 
for convenience to indicate relatively high and in- 
creasing temperature and pressure, but not so high 
that the character of the "elementary" particles 
implied by the choice of L° and h°{T) is destroyed. 
Thus, we assume, for instance, that the pressure is 
not so high as to force an electron from the atomic 
shell, converting the solid into metallic helium (esti- 
mated required pressure 1,000,000 atm), or that the 
temperature is so high as to form a gas mixture of 
ionized particles. This restriction would be of even 
greater importance for diatomic gases, such as hydro- 
gen, which begin to dissociate from molecular to 
atomic particles at relatively low temperatures. 
Presumably, all the properties of the w function 
could be reobtained in such cases by a redefinition of 
L° and h°(T). For example, h°(T) would, then, 
have to be the enthalpy of the gas in the atomic ideal 
gas state and L° the difference between the energy 
of isolated atoms at 0° K and the energy of (for 
hydrogen) the crystal at 0° K and atm. Some of 
these effects, such as dissociation, however, are so 
large that the introduction of a more general w 



function would result in characteristics in the regions 
of interest for cryogenic processes too small to be 
useful. 

From the theoretical point of view this function 
does have the disadvantage that a discontinuity in 
[(dw)/(dP)]r — though not in w — is introduced by 
the liquid-solid transformation for h at 25 atm. 
However, this change in slope of w is slight, being 
given approximately by A[(d/*0/(dP)] T =O.O524j/g- 
atm (from the data of Sw^enson [9]) . It would prob- 
ably be of considerable interest to compare curves 
of w (P Tj T r ) as functions of reduced temperature T T 
and pressure P T for several gases. 

Further insight into the characteristic behavior of 
w(P,T) can be obtained by examination of the 
graphs (at end of paper). From figure 6 which gives 
wasa function of P with T as a parameter, it can be 
seen that all isotherms start with the same value, 
w=\, P = 0; that outside of the liquid II region, all 
isotherms are everywhere decreasing functions of P; 
that the isotherms are continuous except on passing 
through phase changes; that no two isotherms cross; 
that the slopes of all isotherms approach zero with 
increasing pressure; and that, although the total 
percent change in w along an isotherm is greater the 
smaller the temperature, the pressure at which the 
slope of an isotherm, becomes effectively zero in- 
creases with increasing temperature (for the 0° K 
isotherm the slope becomes zero at zero atmosphere, 
the value of w changing discontinuously from one 
to zero through the vapor dome, and then remaining 
zero for all nonzero pressures). 

From figure 7, giving w as a function of T with 
P as a parameter a similar set of statements can be 
made; all isobars start with the same value, w=0 
T=Q;all isobars are everywhere increasing functions 
of T; the isobars are continuous except on passing 
through phase changes ; outside of the liquid II region 
no two isobars cross; the slopes of all isobars outside 
the liquid II region approach zero with increasing 
temperature, and approach zero more rapidly the 
lower the pressure (the slope of the zero atmosphere 
isobar being zero for all nonzero temperatures) ; the 
value of w for all isobars approaches one with in- 
creasing temperature, and approaches more rapidly 
the lower the pressure (the zero atmosphere isobar 
changing discontinuously from zero to one at zero 
atmosphere). The liquid II region sketched on this 
diagram, is somewhat uncertain, due principally to 
lack of time for more extended calculations in an 
area of little interest for liquefier design. 

The limiting behavior of w(p,T) can also be de- 
termined by considering its definition. It is obvious 
that for T-*0, p^0, w^0; for p^0, 2V 0, w->l; 
while for T->0, p=0, or p->0, T=0, w becomes un- 
defined, but lies within the range to 1 . This is the 
region of liquid-vapor coexistence at 0° K and zero 
atmospheres. 

When T increases without limit at constant and 
finite p, 

lim h(p,T) 
Kmw(p,T)= T ^ h0(T) =l. 
r-»«o 
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The second equality of this equation may be 
obtained on the assumption that no matter how high 
the pressure the system will become a "gas" for a 
sufficiently large increase in temperature, with 
enthalpy h approaching the enthalpy of an ideal gas 
h°. This may be seen for a region in which a virial 
equation of state applies, where the higher virial 
type coefficients in the pressure series are sometimes 
taken as decreasing approximately as increasing 
power of 1/T, 

For increasing pressure, however, 



lim w(p,T)= lim , 



i r 

i 5 Jo 



c p (p,T)dT 



=~\ T Vimc p (p,T)dT=0, 

"* Jo p->°° 



since lim c p (2>,T)->0, as can be deduced from the 

p-»co 

fact that the slopes of the isobars with fixed T on a 
temperature entropy diagram, [(i)T)/(dS)] p =(T/c p ), 
approach infinity with increasing pressure (recall, 
however, our restriction on the meaning of P^oo). 

The above combination of well-defined behavior 
and range of magnitude makes the checking and 
cross checking of enthalpy data with w-P and w-T 
diagrams relatively easy. In particular, it has 
indicated many errors and inconsistencies in the 
T-s diagrams of Keeson, and affirmed the consistency 
of the T-s and k-P diagrams of Zelmanov. 

Since the enthalpy of helium can be expressed as 

h(P,T)=h (P)+w(P,T)h°(T), 

the fact that the isotherms on the w-P diagram 
shown in figure (> become approximately horizontal 
at pressures that are not extremely high is of some 
theoretical interest; for this behavior of the iso- 
therms is equivalent to stating that w(P,T) becomes 
practically a function of temperature, alone, even for 
large pressure ranges; thus, since h° is also a function 
of the temperature, alone, the enthalpy can be 
expressed approximately as the sum of a function 
of the pressure h (which might be considered a 
strain energy term) , and a function of the temperature 
w(T)h°(T) (which might be considered a thermal 
energy term). 

3. Liquefier Performance and Heat- 
Exchanger Efficiency 

Before discussing the J-T heat exchanger for 
helium use, it is desirable to review a few general 
relations. Heat exchanger efficiency is defined as the 
ration of actual heat flow to the ideal heat flow, 

^=Cactuai/Qideab where the ideal is that which would 
occur if the product of the total effective heat ex- 
change area and the over-all coefficient of heat 
transfer, AK, is made infinite. However, in obtain- 
ing 77, Q actua i and Qideai must be evaluated for the 
same conditions; that is, the flows, pressures, and 
inlet temperatures of both streams must be held 



constant. Thus, for an arbitrarily defined heat 
exchanger, the efficiency will be an explicit function 
of the flows, and — if the fluid properties are not 
constant — an implicit function of temperatures and 
pressures as well. 

The efficiency of a two-stream, counter-flow, heat 
exchanger can be expressed analytically in three 
different forms, depending on the How conditions. 
If Wi, m 2 are the mass (lows, c pU c p2 , the specific 
heats, Wi = —rhi c pl , W 2 =m 2 c p2 , the heat capacity 
flows in the two streams, and writing, for convenience, 
(}/W) = [(l/W 1 ) + (l/W,d] 9 the expressions for effi- 
ciency are 



-a(1) 
— p \wJ 



^,w>o, 



+ a(T) 



/ _. 1 = — l , w^-w^w- 



(f)(1) >-(!X^ 



1 + 



T )\AK) 



where bars above quantities indicate that some kind 
of integrated average over the heat exchanger must 
be taken. However, obtaining such averages re- 
quires a knowledge of the temperature distributions 
in the streams, and if these were known, there would 
be no problem in the first place. Thus, for heat 
exchangers with variable fluid properties, the above 
formulas can only be used to obtain approximations 
for efficiencies. 

It is our purpose, however, to show how liquefier 
performance depends on an arbitrarily assigned 
efficiency, rather than analyzing how the efficiency 
of a given heat exchanger may be determined. Thus, 
our results will be correct for any exchanger that has 
the assigned efficiency under the assumed flow condi- 
tions. Since the efficiency of a given exchanger 
varies with the flow, an important question is, to 
what extent may a curve of constant efficiency also 
be considered a curve of performance with fixed 
exchanger characteristics? To this we can only give 
a qualitative answer; since the efficiencies of J-T 
heat exchangers are ordinarily very high, even with 
equal heat capacity flows in the streams, the efficiency 
will not be a very strong function of the flows. To 
justify this statement, we use an example, based on 
an approximation of constant and equal fluid proper- 
lies in the streams of a typical counter-flow heat 
exchanger. 

Remembering that the film coefficients are given 
by equations of the form Nu=C Re 08 Pr w , writing 
Kq for the over-all coefficient of transfer for the 
particular case where equal flow exists and <r Q = 
[(AK )/W 2 )] and introducing the liquefied fraction 
a to specify the relation between the mass flows. 
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Wi = — (1 — a) W 2 , the efficiency is given b}^ the 
equation 






1— exp 



2cr a 



-, 0<a<l 



1— a+(l— a) 



or, in the limit as a— >0, as evaluated bv 1'HospitaPs 
Rule, 

i=l+i,a=0. 

n (To 

From these it is easy to see that, for any a, 77 
increases with increasing a ; and, for any <r , 77 
increases with a. However, in the limit, as o- -^oo ? 
the limiting curve 17(a) becomes merely the hori- 
zontal line r)=l. For large, but finite a , 77(a) will 
remain nearly horizontal and nearly straight, so 
that — as mentioned before — 77 is not a very strong 
function of a for J-T heat exchangers. 

To get an idea of the magnitude of the change we 
assume some arbitrary values. Suppose K =50 
w/ft 2 -°K and ,4=25 ft*. Also, let m 2 =0.01 kg/sec, 
^2=5000 j/kg-°K, so that W 2 =50 w/°K, and(7 =25. 
Substitution of these values in the above equations 
shows that the efficiency varies from 96 percent at 
a=0, to 100 percent at a=l. In real liquefiers the 
variation would be even less, since a would ordinarily 
be in the range 0.2 to 0.5. 

4. Calculation of Performance From 
Efficiency 

The liquefaction rate a will now be determined 
from the inlet, high-pressure stream pressure and 
temperature and the heat exchanger efficiency using 
only specific enthalpies from an h—P diagram. In 
making an enthalpy balance, the specific enthalpy h t 
at the inlet of the high-pressure stream is determined 
by the assumed pressure and temperature, and the 
specific enthalpy h g at the inlet of the low-pressure 
stream is assumed to be that of the saturated vapor 
at 1 atm. Thus, only two quantities, the specific 
enthalpy hj above the J-T valve, and the specific 
enthalp}^ h r of the return gas at the outlet of the low- 
pressure stream, are available, and only one of these 
is independent since 

hi—h]=(l — a)(h r —h g ). 

Furthermore, if Ay and h' T are the values of h } and h r 
when 77=100 percent, the efficiency in any other case 
is given by 

h T — h g hi — hj 

h' r — h g h t — h'j 

so that we can write 



and 



hj=hi—r](hi — h f /) 
h r =h g +r)(h r —h' rg ). 



It should be reemphasized here that the efficiency is 
defined in terms of two heat flows (or enthalpy 
changes) under the same conditions of inlet pressures 
and temperatures and the same mass flows. 

Since enthalpies must balance for the whole J-T 
process, as well as just the heat-exchanger, we can 
also write 

h l = ah Iy -\-(l — a)h T 

where h L is the specific enthalp} r of the saturated 
liquid at one atmosphere. Thus a can be expressed by 
either of the equations 

h r — h t h g — h 3 

h r —h L h g —h L 

Substituting for h r and hj their expressions in terms 
of h f r and h'j. and writing L v =h g —h L , these equations 
become 

h t — h g 



77— 



h r r — h g 



h t — h^ hi — h g 



h f r — h g 



Now, a reasonable assumption is that, for 77=100 
percent the temperature difference between streams 
AT becomes zero at some point in the heat ex- 
changer. We should not naively assume that, in the 
most general case, this might not be in the middle of 
the exchanger, since this is exactly what would hap- 
pen if c p for the high-pressure stream were lower 
than that for the low-pressure stream at the higher 
temperatures, and conversely at lower temperatures. 
For helium in the J-T process range, however, the 
reverse is the case, and thus, AT must diminish 
towards the top of the exchanger, and — for suitable 
initial pressure and temperature — diminish again 
towards the bottom of the exchanger. Thus, the 
ideal case will be characterized by having AT=0 at 
either the top or the bottom of the exchanger. If 
we write h' r ' for the value of h' r when A 1=0 at the 
top, and h'/ for the value of h'j when AT=0 at the 
bottom., then — since both h" and h" can be deter- 
mined immediately from an h-P diagram — there will 
be two different possible equations for the lique- 
faction rate a in which all the quantities are known: 



hi—h g 

h' r ' — h g 
a= j ) 

h'r —hg 

hi — h'j' hi — h 



AT-^0 at the top; 



h AT-^0 at the bottom. 



(a) 



(b) a=77 



In case (a) we can determine hj by an enthalpy 
balance and see whether ti^h" ; if not, the as- 
sumption AT=0 at the top must be wrong, and we 
would have to change to case (b). Similarly, in 
case (b) we can determine h' T and check that ti T <Ch" . 
Actually, however, such a check is unnecessary, since 
the two assumptions will result in two curves inter- 
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secting where AT=0 simultaneously at the top and 

hot torn of the exchanger, and it will he immediately 
evident which branch should he used to form the 
composite curve for the liquefaction rate. 

Such composite curves of a as a function of P 
with, rj as a parameter, for helium, are shown in 
figures 9, 10, and 11. The discontinuities in the 
slopes of these curves would be disturbing if one did 
not recall that, even for the case of constant fluid 
properties, the efficiency of an exchanger is expressed 
by two different formulas, depending on whether W 
is positive or negative. It is just the change of the 
total average value of W from, positive to negative 
with the variation of specific heat of helium, which 
produces two independent branches in the curves 
shown. 

The differing way a varies with rj for these two 
branches should be especially noted. For lower 
pressures, relatively large changes in 77 have only a 
small effect on a. If, however, a pressure is chosen 
at or slightly above the "optimum" pressure, a will 
not only decrease rapidly with increasing pressure, 
but at a given pressure, will drop drastically with 
slight reduet ions in beat-exchanger efficiency. There- 
fore, it is desirable to be somewhat conservative in 



the design pressure for a liqucfier in order to allow 
for possibly too optimistic estimates of heat-ex- 
changer performance. 
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Figure 1. Specific enthalpy of helium. 
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Figure 2. Specific enthalpy of helium at 0° K, h (p) versus 
pressure. 
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[(d/i)/(dp)]r versus pressure of helium below 4° K. 

from density and expansivity data from Keesom, helium. 
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Figure 3. Specific volume of helium at 0° K. 
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Figure 4. Compressibility of helium at 0° K. 
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Figure 6. Helium w-p diagram. 
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Figure 7. Helium w-T diagram. 




Figure 8. Yield a as a function of pressure and temperature at the J-T valve. 




Figure 9. Yield a with J-T heat exchanger efficiencies rj equals 100 percent as a 
function of HP inlet temperature and pressure. 
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Figure 10. Yield a. with J -T heat exchanger efficiency t\ equals 95 percent as a 
function of HP inlet pressure and temperature. 
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Figure 11 Yield a with J-T heat exchanger efficiency y equals 90 percent 
as a function of HP inlet pressure and temperature. 
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Figure 12. Yield a. for several heat exchanger efficiencies as a function oj inlet 
temperature and pressure. 



Boulder, Colo., August 30, 1957. 
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